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The analysis of the terminating bands has been performed in the relativistic mean field framework. 
It was shown that nuclear magnetism provides an additional binding to the energies of the specific 
configuration and this additional binding increases with spin and has its maximum exactly at the 
terminating state. This suggests that the terminating states can be an interesting probe of the 
time-odd mean fields provided that other effects can be reliably isolated. Unfortunately, a reliable 
isolation of these effects is not that simple: many terms of the density functional theories contribute 
into the energies of the terminating states and the deficiencies in the description of those terms 
affect the result. The recent suggestion [1] that the relative energies of the terminating states in the 
N =f= Z, j4 ~ 44 mass region given by AiJ provide unique and reliable constraints on time-odd mean 
fields and the strength of spin-orbit interaction in density functional theories has been reanalyzed. 
The current investigation shows that the AE value is affected also by the relative placement of 
the states with different orbital angular momentum I, namely, the placement of the d {I = 2) and 
/ (Z = 3) states. This indicates the dependence of the AE value on the properties of the central 
potential. 

PACS numbers: PACS: 



I. INTRODUCTION 

The density functional theory (DFT) in its non- 
relativistic [2] and relativistic [3] realizations is a stan- 
dard tool of modern nuclear structure studies. However, 
providing global description of atomic nuclei, it still suf- 
fers from the fact that many channels of effective inter- 
action are not uniquely defined: this is a reason for a 
large variety of the DFT parametrizations, the quality 
of many of which is poorly known. The spin-orbit in- 
teraction and the time-odd mean fields are of particu- 
lar interest in this context, since there are considerable 
variations for these quantities (see, for example, Refs. [4- 
6]). The spin-orbit interaction plays a crucial role in the 
definition of the shell structure of nuclei, and, thus its 
accurate description is required so that theoretical tools 
have predictive power for nuclei beyond known regions. 
The time-odd mean fields (or nuclear magnetism (NM) in 
the language of the relativistic mean field (RMF) theory 
[7]) contribute to the single-particle Hamiltonian only in 
situations where the intrinsic time-reversal symmetry is 
broken and Kramers degeneracy of time-reversal counter- 
parts of the single-particle levels is removed. The rotat- 
ing nuclei and odd and odd-odd mass nuclei are typical 
examples of such situations, see Refs. [2, 3] and references 
quoted therein. 

It was recently suggested in Ref. [1] that the set of 
terminating states in the N Z, A ~ 44 mass re- 
gion provides unique and reliable constraints on time- odd 
mean fields and the strength of spin- orbit interaction in 
Skyrme density functional theory (SDFT), see also Refs. 
[8, 9]. Later this procedure (called as 'TS-method' in this 
manuscript, where 'TS' refers to 'terminating states') has 
been used in the analysis of terminating states in this 
mass region within the RMF theory [10] which is one 
of the versions of covariant density functional (CDFT) 



theory [3]. 

The authors of Refs. [1, 8-10] claim that the TS- 
mcthod is free from the drawbacks of standard methods 
of defining spin-orbit interaction based on measuring the 
single-particle energies of the spin-orbit partner orbitals 
in spherical nuclei. As a consequence, it is stated that it 
allows to define very accurately both isoscalar and isovec- 
tor channels of spin-orbit interaction [1, 10]; the feature 
which was impossible in the previously existing methods. 

The conclusions obtained within the TS-mcthod are 
drastically different from the ones previously obtained in 
the SDFT and RMF frameworks. For example, based on 
the comparison of the calculated and experimental ener- 
gies of spin-orbit partner orbitals, it was shown in Ref. 
[5] that the experimental spin-orbit splittings are bet- 
ter reproduced in the RMF approach than in the SDFT 
(see Fig. 2 in Ref. [5]). On the contrary, the results ob- 
tained in Refs. [1, 10] within the TS-method show that 
the SDFT provides better description of spin-orbit split- 
tings than the RMF: it was suggested in Ref. [1] that 
only 5% reduction in isoscalar Is strength is needed in 
the SDFT approach in order to reproduce experimental 
data. Considering the confiict of these results and the 
importance of the spin-orbit interaction in nuclei it is 
necessary to understand to which extent the basis of the 
suggested TS-method is sound and justified. The goal 
of the present manuscript is the study of the properties 
of terminating bands and their terminating states in the 
RMF framework. In particular, the question of whether 
all DFT contributions have been correctly accounted in 
the realization of the TS-method in Refs. [1, 10] is ad- 
dressed is the current manuscript. 

The manuscript is organized as follows. Time-odd 
mean fields in terminating bands are studied in Sect. 
II. The basis of the TS-method, its realization in self- 
consistent DFT and in the Nilsson potential are discussed 
in Sect. III. Sect. IV analyses the contributions of differ- 
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FIG. 1: (Color online) Calculated total angular momentum 
(panel (a)), total binding energy (panel (b)), kinematic and 
dynamic moments of inertia (panels (c) and (d)) as a func- 
tion of rotational frequency i}x in the ground state band of 
^"Ne. The results obtained with and without nuclear mag- 
netism are presented. Note that the band termination takes 
place at rotational frequency above which the subsequent in- 
crease of rotational frequency does not modify neither total 
angular momentum nor total binding energy. This point also 
corresponds to 7 = 60° (Fig. 3c). 



ent DFT terms into the relative energies of terminating 
states in the A ~ 44 mass region. The discussion of the 
energy scale, its connection to the effective mass of the 
nucleon and their impact on the relative energies of ter- 
minating states is presented in Sect. V. Finally, Sect. VI 
summarizes main conclusions. 



II. TIME-ODD MEAN FIELDS IN 
TERMINATING BANDS: TEST CASE OF ^onE 

Previous DFT investigations of the modifications of the 
moments of inertia [4, 11-13] and single-particle proper- 
ties [6] in rotating nuclei caused by the time-odd mean 
fields (nuclear magnetism) were restricted to the superde- 
formed (SD) bands. However, these bands are far away 
from the termination and are characterized by a rela- 
tively stable deformation. In order to understand how 
NM affects the properties of the terminating bands, the 
ground state configuration in ^°Ne has been studied. 
This band is a classical example of band termination [14]. 
It has the TT{d5/2)it^{d5/2)'i configuration relative to the 
^^O core with maximum spin Imax = 8"*". The selection 
of this configuration has been guided by its simplicity, 
which allows us to understand the role of time-odd mean 
fields in terminating bands in greater details. Although 
the terminating bands were observed also in heavier nu- 
clei, it is difficult to trace them from low spin up to band 
termination in the self-consistent approaches [3] without 
special techniques such as used in the cranked Nilsson- 
Strutinsky (CNS) approach [14]. 



The investigation of NM in ^"Ne has been carried out 
within the framework of the cranked relativistic mean 
field (CRMF) theory [12, 15] following the formalism of 
Ref. [6], where similar study has been performed for the 
yrast SD band in ^^^Dy. In the CRMF calculations of 
this manuscript, all fermionic and bosonic states belong- 
ing to the shells up to Np = 12 and Nb = 16 are taken 
into account in the diagonalization of the Dirac equa- 
tion and the matrix inversion of the Klein- Gordon equa- 
tion, respectively. The detailed investigation indicates 
that this truncation scheme provides very good numeri- 
cal accuracy. The NLl [16] parametrization of the RMF 
Lagrangian is employed in the study of ^'^Ne, while the 
studies of terminating states in the A ^ 44 mass region 
(Sects. Ill and IV) are performed mostly with the NLl 
and NL3 [17] parametrizations. The pairing is neglected 
in calculations. 

Fig. 1 shows the total angular momentum, total bind- 
ing energy E and kinematic (J^^^) and dynamic (J*-^') 
moments of inertia of the ground state configuration in 
^•^Ne as a function of rotational frequency obtained in the 
calculations with and without NM (the later will be fur- 
ther denoted as WNM). The band crossing caused by the 
interaction of the r = +i signatures of the [220] 1/2 and 
[2 11] 3/2 orbitals both in the proton and neutron subsys- 
tems takes place at lower frequency in the calculations 
with NM; this is in line with previous finding that the 
NM shifts the band crossing frequencies [18]. The NM 
also increases the moments of inertia before band crossing 
(Fig. lc,d); similar effect has been seen before in the SD 
bands (see Ref. [6] and references therein). It also leads 
to a faster alignment of angular momentum with rota- 
tional frequency (Fig. la); full alignment at Imax = 8"'' 
corresponding to band termination takes place at lower 
frequency in the calculations with NM. 

In the context of study of terminating states two re- 
sults are important. First, at the band termination the 
NM does not modify neither total angular momentum 
(Fig. la) nor the expectation values of the single-particle 
angular momenta < >,; (Fig. 2). At lower frequency, 
the impact of NM on < jx >i is similar to the one previ- 
ously studied in the SD band of ^^^Dy [6], and, thus, it 
will not be discussed in detail. However, one should men- 
tion that when analyzing the impact of NM on < >i, 
the region of band crossing and the region close to the 
band termination have to be excluded from considera- 
tion because considerable differences in the deformations 
of the NM and WNM solutions at given frequency distort 
their comparison. Second, the NM provides an additional 
binding to the energies of the specific configuration and 
this additional binding increases with spin and has its 
maximum exactly at the terminating state (Fig. lb and 
3d)). This suggests that the terminating states can be 
an interesting probe of the time-odd mean fields provided 
that other effects can be reliably isolated. 

When the results of the NM and WNM calculations 
are compared as a function of total angular momentum, 
one can see that the quadrupole deformation (32 (Fig. 
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FIG. 2: (Color online) The expectation values of single- 
particle angular momenta < jx >i of the neutron orbitals 
occupied at low frequency in the ground state configuration 
of 2°Ne given along the deformation path of this configura- 
tion. The single-particle orbitals are labeled by means of 
the asymptotic quantum numbers [NnzA]Q (Nilsson quantum 
numbers) of their dominant component of the wave function 
at fia; — 0.0 MeV. The orbitals with signature r = +i and 
r = —i are shown in right and left panels, respectively. The 
results of the calculations with and without NM are shown by 
solid and dashed lines, respectively. The region of the band 
crossing is located between the dashed lines. 



3a), mass hexadecapole moment Q40 (Fig. 3b), and 7- 
deformation (Fig. 3c) are almost the same in both calcu- 
lations. The only difference is seen in the total binding 
energies (Fig. 3d) , where the NM solution is more bound 
than the WNM solution. These results give a hint why 
the cranked models based on the phenomenological po- 
tentials like Woods-Saxon or Nilsson, which do not in- 
clude time-odd mean fields [4], are so successful in the 
description of experimental data. When considered as 
a function of spin the deformation properties of the ro- 
tating system are only weakly affected by the time-odd 
mean fields, and the proper renormalization of the mo- 
ments of inertia [14] takes care of the E versus angular 
momentum curve. 



0.6 
^0.5 

c 0.4 

I 0-2 
Q 

0.1 
60 



E 
o 

% 20 



■ WNM 

■ NM 



(C) 



(b) 



(d) 



150 
100 



a 



-50 
-147 



4 6 8 2 4 6 
Angular momentum (h) 



-149 
-151 
-153 
-155 



FIG. 3: Calculated quadrupole deformation P2 (panel (a)), 
mass hexadecapole moment Q40 (panel (b)), 7-deformation 
(panel (c)), and total binding energy E (panel (d)) as a func- 
tion of angular momentum. The results obtained with (NM) 
and without (WNM) nuclear magnetism are presented. 



III. THE TS-METHOD 

The TS-method suggested in Ref. [1] employs the ter- 
minating states in the A ~ 44 mass region with the 
proton /"/^^ ^^'^ /7/2 structure, but, in general, 
according to Ref. [9] can be employed to terminating 
states in any mass region. The difference AE^'^p = 
-^(^3/2 f7/2^) ~ ^(/7/2) between the excitation energies 
^(•^3/2 /"/V) ^^"-^ ^(17/2) "-"^ terminating states with 
the structure d~J^^ /y/t^ -^7/2 dominated by the size 
of the magic gap 20 which is surrounded by the ^3/2 and 
/7/2 spherical subshells (see Fig. 4). 



A. The spin-orbit splittings in the TS-method 

The principal difference between the standard and the 
TS-methods of defining the strength of spin-orbit inter- 
action is schematically shown in Fig. 4. The standard 
method requires that both partners of spin-orbit Is dou- 
blet with j = I ± J are observed in experiment since the 
spin-orbit splitting AE^'~' is related to the strength of 
spin-orbit interaction. This severely restricts the possi- 
bilities to study spin-orbit interaction since both part- 
ners should be located in the vicinity of the Fermi level 
to be observed: this condition is very difficult to sat- 
isfy for high-j orbitals since they are characterized by 
large spin-orbit splittings, see, for example, Refs. [5, 19]. 
On the contrary, the TS-method employs the terminat- 
ing states based on the particle-hole excitations involving 
the single-particle states with j = I — j and / = i' + j 
which emerge from different A^-shells and are character- 
ized by the energy splitting AE'^^ (Fig. 4). Since these 
states are located in the vicinity of the Fermi level, the 
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TS-method provides also information on the spin-orbit 
interaction of high-j orbitals according to Ref. [1]. 

It is necessary to recognize that both methods of defin- 
ing the spin-orbit interaction are not free from important 
drawbacks. The experimental single-particle states in 
spherical nuclei used in the standard method are strongly 
affected by the couplings with vibrations in many cases 
[20]. On the other hand, the AE^^ value used in the TS- 
method depends not only on the spin-orbit splitting but 
also on how well the positions of the single-particle states 
with different orbital momenta I and V (Fig. 4) are de- 
scribed in the DFT calculations. The later fact has been 
neglected in Ref. [1] using an analogy with the Nilsson 
potential, the validity of which is questioned below. 



• the polarization effects (both in time-even and 
time-odd channels) on going from the fy^ to the 

'^3/2 /"/V terminating state. 

For simplicity of discussion, they will occasionally be 
called as 'ingredients of Ai?'^'^'. 

can also be split into the terms which depend 
on time-even (TE) and time-odd (TO) mean fields 
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The minus sign in front of ^Ej,q reflects the fact that 
NM always decreases the size of /S.E. 
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FIG. 4: Schematic comparison of the standard and TS- 
methods of defining the strength of spin-orbit interaction. In 
the left and right panels, the single-particle spectra without 
(on left) and with (on right) spin-orbit interaction are com- 
pared. 



B. The TS-method in self-consistent approaches 



In the self-consistent calculations, the AE^^ = 
^('^3/2 /"/t^) ~ ^(^7/2) quantity is defined as the differ- 



'7/2 



ence of the binding energies of the corresponding termi- 
nating states. Without going into the details of specific 
DFT (nonrelativistic or relativistic), one can conclude 
that /S.E^^ depends on 

• the energy scale of the single-particle spectra which 
is related to the effective mass m*{kp)/m of the 
nucleon at the Fermi surface, 

• the spin-orbit interaction, 

• the relative placement of states with different an- 
gular momentum Z, 



The ingredients of AE^^ depend in a complicated way 
on different terms of the DFT with at least one term 
contributing into each of four first ingredients of AE^'-^ 
within the nonrelativistic SDFT. In the RMF theory, the 
spin-orbit interaction is defined in a natural way with- 
out additional coupling constant [3] . The time-odd mean 
fields related to NM are defined through the Lorcntz co- 
variance [3] and also do not require an additional coupling 
constant. However, both these terms depend in an indi- 
rect way on the coupling constants of other terms of the 
RMF Lagrangian. 

Considering the uncertainties of the description of the 
'ingredients of AE^^^ in the DFT, it is not obvious that 
simple fit (to experimental AE'^^p values) of the coupling 
constants of the DFT terms related to a pair of ingredi- 
ents of AE^'-^ (such as time-odd mean fields and spin- 
orbit interaction as in Ref. [I]; the later treated pertur- 
batively) will allow to define these constants in a unique 
way. This is especially true considering that physical ob- 
servables depend on many (or sometimes all) coupling 
constants simultaneously within the DFT, and the effect 
of varying one or two coupling constants may be either 
enhanced or cancelled by a variation of others. Strictly 
speaking, the quality of such perturbative fits involving 
only one or two terms of the DFT is not known until 
the results of global fit including all the DFT terms are 
available. For example, it was shown in Ref. [21] that 
perturbative studies of tensor terms allow only very lim- 
ited conclusions. 



2. Polarization effects 

Fig. 5 illustrates the polarization effects present in the 
CRMF calculations of terminating states. Since fully 
stretched states with spin Imax are reasonably well de- 
scribed by a single Slater determinant [23], the compar- 
ison with experiment is performed without angular mo- 
mentum restoration as in all DFT studies of the termi- 



nating /„ 

the time-odd mean fields (see Sects. II and IV B), pie, Refs. 



ax states in this mass region, see, for exam- 
[1, 10]. At spherical shape, the f'j/2 — ^3/2 
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FIG. 5: (Color online) (right part) Single-particle energies 
at spherical shape in ''^Ti obtained in the RMF calcula- 
tions with the NLl parametrization and the Nilsson po- 
tential with standard set of parameters [22] (columns indi- 
cated as "RMF-spher" and "Nilsson"). In order to facili- 
tate the comparison, the RMF states are given at the cal- 
culated energies, while the energies of the states obtained in 
the Nilsson potential are shifted in such a way that the av- 
erage energy of the I/7/2 and 1^3/2 states is the same in 
both calculations, (left part) The magnitudes AE of the 
/7/2 — ^3/2 splittings as extracted from the energies of ter- 
minating states {AE = E{d^J^fj^^) — E{fy2)) are shown 
by arrows in columns "exp" (experiment), "CRMF-NM" and 
"CRMF-WNM" (CRMF calculations with and without NM). 
The middle points of arrows are located at the average en- 
ergy of the 1/7/2 and ld3/2 RMF states. The terminating 
states are deformed in the CRMF calculations. Note that the 
value of A_E = 5.1 MeV (not shown in figure) obtained in the 
cranked Nilsson-Strutinsky calculations compares favorably 
with experiment [AE"'''' = 5.51 MeV). 



difference between experiment and calculations consid- 
erably (by 1.14 MeV) (column "CRMF-NM" in Fig. 5). 
Note that mass and charge quadrupole and hexadecapole 
moments change only by ^ 10~^% on going from the 
CRMF-WNM to CRMF-NM solutions. Thus, the de- 
formation differences between these two solutions are al- 
most non-existant and can be neglected. Based on the 
consideration of polarization effects, the AE^*-^ can be 
approximated as 



AE^'^ « AE^P^ + AE^'-'^^^P"'' 



AE^o 



(2) 



Considering that the terminating states of interest are 
close to spherical (Fig. 6), this approximation which cor- 
responds to a perturbative treatment of the deformation 
polarization effects should be quite reasonable. This ap- 
proximation also allows to use the results of spherical 
RMF calculations in subsequent analysis of AE (Sect. 
IV). 
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FIG. 6: (Color online) Quadrupole deformations of the ter- 
minating states obtained in the CRMF calculations with the 
NLl parametrization. For each nucleus, the deformation po- 
larization energies Ai?'*'^-'"''"' are shown between the cal- 
culated deformation points. Three nuclei with the largest 
^]7;'ief-poi yg^jygg g^j-g indlcated by arrows. 



energy gap {AE'^p'^) in the single-particle spectra consid- 
erably exceeds AE'''=p. The AE'^p'' value is a very good 
approximation to the results of the spherical RMF cal- 
culations without NM in which the gap between these 
states is defined as the difference of the binding ener- 
gies of the c?^y2 /"/V ^^'^ f7/2 states: the difference be- 
tween two results for all nuclei under study does not ex- 
ceed 40 keV. When deformation polarization effects de- 
noted as AE'^''f~P°^ are taken into account (the "CRMF- 
WNM" column in Fig. 5), this gap becomes even larger 
and reaches AE = 7.57 MeV exceeding by 37% the ex- 
perimental value. The inclusion of NM decreases the 



C. The Nilsson potential analogy of Ref. [1] 

In order to overcome the problems discussed in Sect. 
IIIB 1, the authors of Ref. [1] use the analogy with the 
simple form of the Nilsson potential [24] for the AE^'^ 
part. This form of the Nilsson potential is similar to 
the one given below (Eq. (5)), but with the parameters 
K and fi independent on principal quantum number TV. 
However, no proof (analogy is not a proof) is provided 
whether such a transition from self-consistent DFT to the 
Nilsson potential is valid and whether the dependence 
of the energies of the single-particle states on quantum 
numbers N,j,l, and s is the same in the self-consistent 
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DFT and in the Nilsson potential. Fig. 5 clearly shows 
that the later is not a case. 

In the simple form of the Nilsson potential, the magni- 
tude of the /7/2 — (^3/2 splitting related to the magic gap 

20 is given by 

A£;^*' =fejo(l-6K-2KAt). (3) 

Thus, it depends on three major factors: (i) the en- 
ergy scale of the single-particle potential characterized 
by fkJo, (ii) the flat-bottom and surface properties en- 
tering through the orbit-orbit term ^, and (iii) the 
strength of the spin-orbit term k. Then, the authors of 
Ref. [1] using the fact that in light nuclei the Nilsson po- 
tential resembles the pure harmonic oscillator potential, 
which leads to /z ~ 0, conclude that the magnitude of the 
h/2 - d3/2 splitting is given by 

AE^'^ = /iWo(l - 6k). (4) 

Thus, in this approximation, the A£^*' splitting is de- 
fined only by the energy scale TujOq and the strength of 
the Zs-potcntial. Arguing that the energy scale is rather 
well constrained by the data not only for the Nilsson but 
also for the self-consistent approaches, the authors of Ref. 
[1] conclude that the /7/2 — c?3/2 splitting is directly re- 
lated to the strength of the Zs-potential. Or alternatively, 
this approximation corresponds to the situation when the 
^j^Niis (jQgg jjQ^ depend on orbital motion of the nucle- 
ons. 



D. An alternative form of the Nilsson potential 

One question of paramount importance we have to ask 
is whether the simple form of the Nilsson potential given 
in Ref. [24] and used in Ref. [1] is unique and how well 
it describes the experimental data. It turns out that 
the modern versions of the Nilsson potential employ the 
parameters k and ji which are dependent on the main os- 
cillator quantum number A'' and on nucleon type (proton 
or neutron), thus facilitating the study of wide range of 
nuclei with the same set of single-particle parameters and 
with comparable accuracy [22, 25, 26]. Since the Nilsson 
potential is phenomenological in nature, this procedure is 
well justified. The accuracy of the description of different 
physical quantities such as, for example, rotational prop- 
erties and relative energies of different single-particle con- 
figurations [14, 22, 25, 27] is considerably improved when 
different values of k and /x are used for different A''-shells. 
Although some variations between the parametrizations 
exist, this approach is used in almost all parametriza- 
tions of the Nilsson potential developed from the middle 
of the SOtics of the last century [22, 27, 28]. The studies 
employing this description of the Nilsson potential are 
abundant and provide systematic information on the ac- 
curacy of the description of physical observables in differ- 
ent mass regions. Even more sophisticated dependence of 



TABLE I: The standard parametrization of the Nilsson po- 
tential from Ref. [22, 25]. Only the parameters for the N = 2 
and 3 shells are shown. 
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0.105 


0.00 


0.105 


0.00 
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0.090 


0.30 


0.090 


0.25 



the K and ji parameters on the principal quantum num- 
ber A'' and orbital angular momentum I is introduced in 
Ref. [29] and employed in a number of studies, see, for 
example, Refs. [30, 31]. 

1. Terminating states in the A^^j 4A mass region 

In the Nilsson potential with kjv and /ijv parameters 
dependent on the principal oscillator quantum number 
N [22, 25] 

jjNil-N-dep _ ^ 
= hlUo{N-KN[2ls + ^lN{f-<f>N)]}, (5) 

the magnitude of the — sphtting related to the 
energy difference of the d'^J^ /vyt^ and fy^ terminating 
states is given by 

= hUo{l - 3[k2 + Ks] - 3K3/i3 - K2M2)- 

The superscript 'Nil — N — dep' is used to indicate the 
dependence of these expressions on the main oscillator 
quantum number N. 

The kn and parameters of the so-called standard 
parametrization of the Nilsson potential are shown for 
the shells of interest in Table I. For the protons, the 
value AE^''^-^-'^''P = 0.415S.cjo is obtained in the cal- 
culations employing the kn parameters from the Ta- 
ble I but assuming the hn = as it was done in the 
derivation of Eq. (4) . This value can be compared with 
the AE^''--^-'^''P = 0.334fujo value obtained with the 
kn,I^n parameters from the Table I. One can see that 
these two values differ by approximately 25% and this dif- 
ference is solely attributed to the non-zero value of the 113 
parameter. Considering that AE'^'^p ~ 5.5 MeV (see Fig. 
7), 25% difference correspond to 1.4 MeV; this difference 
definetely cannot be ignored when experimental data is 
compared with experiment. 

2. Concluding remarks 

Even for terminating states in the A 4:4: mass re- 
gion one cannot ignore the dependence of the energies 
of the (N,l) and (N',l') states, from which the j = I — j 
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and j' = I' + J states used in the TS-method emerge (see 
Fig. 4), on the orbital angular momentum. This depen- 
dence enters through the iJ,]sr{P— < P >n)\ term of the 
Nilsson potential (Eq. (5)). This is contrary to the ap- 
proximation made in the derivation of Eq. (4) which has a 
consequence that the energy difference A£^^*'* depends 
only on the energy scale hwf) and the strength of the 
spin-orbit term k. The dependence of AE'^*''^ on the or- 
bital angular momentum in the case of terminating states 
involving single-particle states from higher iV-shells has 
been recognized in Ref. [9]. 



IV. TERMINATING STATES IN THE ^ ~ 44 
MASS REGION: WHAT WE CAN LEARN FROM 
THE COMPARISON WITH EXPERIMENT 
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FIG. 7: (Color online) The experimental and calculated mag- 
nitudes Ai5 of the fi/2 — d-1,/2 splittings as extracted from the 
energies of terminating states (Ai? = E{d~j^ /"/t^) — ^(/7/2))- 
The results of the CRMF calculations are shown for the NLl 
and NL3 parametrizations of the RMF Lagrangian. The ex- 
perimental data corrected for the deformation polarization 
effects (as obtained in the NLl parametrization) is shown by 
open circles. Note that the experimental and deformation 
polarization corrected values of Ai5 coincide in the case of 

47y^ 



Fig. 7 compares the results of calculations with exper- 
iment. The same data set as in Refs. [1, 10] is used in 
this comparison, but it is shown as a function of n, where 
n stands for the number of the f'j/2 protons in the fy^ 
terminating state. In addition, this figure compares the 
absolute values and not the differences between experi- 
mental results and calculations as it was done in Refs. 
[1, 10]: the differences normalized to ''^Ca are compared 
in Fig. 8a. Since the n value is the same for each iso- 
tope chain (n = for the Ca isotopes, n = 1 for the 
Cs isotopes, n = 2 for the Ti isotopes, and n = 3 for 
the V isotopes), the isospin dependence of AE is clearly 



visible. Different isotope chains show different isospin 
dependencies and they are well reproduced in the cal- 
culations (Fig. 7 and Fig. 8a). On the other hand, the 
calculations overestimate the absolute value of AE (Fig. 
7), and the difference between the calculated and exper- 
imental AE values show pronounced n-dependence (Fig. 
8a). 

One source of the discrepancy between theory and ex- 
periment is related to the impact of effective mass of the 
nucleon on the single-particle spectra (see Sect. V): in 
general, it should lead to an overestimate of experimen- 
tal AE in the calculations. The other sources of these 
discrepancies are analyzed in detail in this Section. 
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FIG. 8: (Color online) (a) The difference (A£;*''-A.B*=^f )„<,rm 
between the calculated and experimental values of AE (based 
on the results of Fig. 7). This difference is normalized to 
zero for **Ca. (b) The calculated difference ETo{d^j2 /"/V)" 
^To{.f 7/2) shown as a function of n for the indicated RMF 
parametrizations (based on the results of Fig. 9). 



A. Deformation polarization effects 



The deformation polarization effects discussed in Sect. 
1IIB2 are characterized by the AE'^'^^~'p°^ energies. The 
sign of AE'^'^^^P"' depends on relative deformations of 
the f^,2 ^3/2/7/2^ terminating states (Fig. 6). It 



is positive (negative) when the /32-deformation of the 
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"^3/2/7/2 state is larger (smaller) than the one of the 
/^/2 state. The AE'^^f-P"^ values almost do not depend 
on the parametrization of the RMF Lagrangian: the dif- 
ference in their values is below 10 keV if the results of 
the NLl and NL3 parametrizations are compared. If the 
experimental data are corrected for these deformation 
polarization effects, then smooth trend (the curve 'exp. 
(def-cor)' in Fig. 7) as a function of n emerges. The 
AE value along this curve decreases by ^ 0.25 MeV on 
going from rt = to n = 3. Assuming that these effects 
are reasonably well described in the calculations, one can 
conclude that the nucleus-dependent fluctuations in ex- 
perimental value of AE (the curve 'exp.' in Fig. 7) are 
due to deformation polarization effects. 



B. Nuclear magnetism (time-odd mean fields) in 
the terminating states of the A ^ 44 mass region 

Fig. 9 shows the additional bindings Eroistate) to the 
energies of terminating states due to NM. This quantity 
increases with the increase of the n- value for the ^-^.d 

•^3/2 /-T/t^ terminating states. The increase of Eto with 
isospin within specific isotope chain is associated with 
the increase of the number of the occupied neutron /y /2 
states and corresponding increase in spin. The increase 
of the values of Eto correlates with the increase of the 
spin of the terminating states: for example, the fy^ and 

^3/2 terminating states have Imax = 6^ and Imax = 



= ^ and /„ 



35 + • 

= f m 



V, 



11^ in Ca and Imax 
respectively. 

The results of the calculations confirm previous con- 
clusion obtained in ^°Ne (Sect. II) that the additional 
binding due to NM is considerably enhanced in the ter- 
minating states. At no rotation, the additional binding 
due to NM to the energies of the single-particle config- 
urations in odd-mass nuclei is in average around ~ 100 
keV and seldom reaches 200 keV in the mass region of 
interest [32]. This is much smaller than the additional 
binding observed in the terminating states in which it 
reaches 4 MeV for n = 3 (Fig. 9). 

Because of their magnitude, the Eto values in ter- 
minating states are also a good measure of how well the 
time-odd mean fields are defined in the specific version of 
DFT. The Eto values obtained with different frequently 
used non-linear parametrizations of the RMF Lagrangian 
such as NLl [16], NL3 [17], NLSH [33], NLRAl [34] and 
NLZ [35] are shown in Fig. 9. With increasing Eto and 
n, the absolute variations in the Eto values calculated 
with different RMF parametrizations increase. However, 
they are still within 15% of the absolute value of Eto- 
This result suggests that within the non-linear versions 
of the RMF Lagrangian NM is defined with similar ac- 
curacy. 

This value can be used to estimate the uncertainty in 
the definition of the moments of inertia in the CRMF 
calculations due to the uncertainty in NM. Dependent 
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FIG. 9: (Color online) Additional bindings Eto (state) (in 
absolute value) to the energies of terminating states due to 
NM shown for the terminating states of interest. The results 
are shown for the indicated parametrizations of the RMF La- 
grangian as a function of n. 



on the nuclear system, the NM contribution to the to- 
tal kinematic moment of inertia is approximately 10-25% 
[12, 32]. Thus, the uncertainty of the definition of the 
absolute value of the total kinematic moments of inertia 
due to the uncertainty in the definition of NM is modest 
being in range of 1.5-3.75%. 

It follows from Fig. 8 that the portion of the n- 
independent part of the discrepancy between experimen- 
tal and calculated AE values may be related to the un- 
certainties in NM since the difference EToid^/2 /t/V) ~ 
ETo{fy2) somewhat (within « 200 keV) depends on the 
RMF parametrization. The contribution of NM into the 
n-dependent part of this discrepancy is discussed in Sect. 
IV C 3. 



C. The dependence of Ai5 on orbital angular 
momentum and spin-orbit interaction 

1. The impact of density modifications on the 
single-particle properties 

It is well known fact that the modifications of the cen- 
tral nucleonic potential and its surface properties affect 
the single-particle states with different angular momen- 
tum I in a different way (see Refs. [20, 36] and references 
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therein). They also alter the spin-orbit potential and 
lead to the changes in the spin-orbit splittings. In order 
to check how big this effect is in the nuclei under study, 
the proton density distributions and single-particle spec- 
tra at spherical shape are compared in Fig. 10 for the 
'^^Ca. and ^^V nuclei. These nuclei represent the lower 
and upper mass ends of the data set under investigation. 
The configuration of V has 2 additional /y /2 neutrons 
and 3 additional /y /2 protons as compared with the con- 
figuration of ^^Ca. The filling of these high-j orbitals 
increases the density near the surface (Fig. 10a). These 
modifications of the density change the central and spin- 
orbit nucleonic potentials (in a similar fashion as it was 
discussed in Refs. [37, 38]) leading to the modifications 
of the single-particle spectra (Fig. 10b). 

On going from *^Ca to ^^V, the spin-orbit splitting in 
the dr^ji — d-iji doublet decreases by 0.12 MeV (from 6.72 
MeV to 6.60 MeV), while the one in the J-; ji — dou- 
blet increases by 0.57 MeV (from 7.0 MeV to 7.57 MeV). 
If the modifications in the single-particle spectra would 
be restricted only to the spin-orbit splittings and their 
modifications would evenly be redistributed between the 
i = 1 + 1/2 and j = 1 — 1/2 members of the spin-orbit 
doublet, this would decrease the I/7/2 — 1^3/2 splitting 
by 0.22 MeV. 

However, the calculations show that the fjf2 ^ '^3/2 
splitting is increased by 0.34 MeV (Figs. 10). Assum- 
ing that the energy scale does not change on going from 
"'■^Ca to ''^V, this can only be explained by the change 
of the relative positions of the d {I = 2) and f {I = 3) 
states from which the 6^3/2 and fj^2 states emerge. Un- 
fortunately, there is no straightforward way in the RMF 
calculations to get an access to the {N, I) states (in sense 
of Fig. 4). Thus, in order to illustrate the dependence on 
I, the centroid energy (denoted as "centr(state)" in Fig. 
10b) and defined as an average energy of the members of 
spin-orbit doublet is used. Fig. 10b shows that the energy 
gap between the centroids of the d and / spin-orbit dou- 
blets increases by 0.55 MeV on going from "^^Ca to "^^V. 
As a consequence of this increase and the above discussed 
changes in the spin-orbit splittings, the I/7/2 — 1^3/2 
splitting increases by 0.34 MeV. This value represents 
more than half of the increase of (Ai?*'' — A_E*^^^')„orm 
on going from *^Ga to ^^V (Fig. 8a). 



Relative placements of the states with different angular 
momentum I 
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FIG. 10: (Color online) (a) Proton density distribution in 
the ''^Ca and nuclei as obtained for ground state in the 
spherical RMF calculations with the NLl parametrization of 
the Lagrangian. (b) Corresponding energies of the single- 
particle states. The states in ''^Ca are shown at the calculated 
energies, while all states in ''^V are shifted by constant value 
in such a way that the average energy of the 1^3/2 and I/7/2 
states is the same as in *^Ca. 



The fact that the relative placement of states with dif- 
ferent orbital angular momentum I (especially, of high- 
l states) is not well reproduced in non-relativistic and 
relativistic mean field models is well known, see Refs. 
[20, 21, 36]. The origin of this problem is connected 
with the surface profile of the mean field and kinetic 
terms. Microscopic considerations indicate that the effec- 
tive mass of the nucleon has a pronounced surface profile 
which is insufficiently parametrized in the present mean 



field models [20]. In Refs. [1, 10], this fact has been ig- 
nored and no proof has been provided that the placement 
of the d and / states, from which the ^3/2 and /7/2 states, 
emerge is correct. 

It turns out that the difference in absolute value of AE 
obtained in the NLl and NL3 parametrizations (Fig. 8) is 
well explained by the differences in the relative energies 
of the d and / states in these parametrizations. The 
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energy gap between the ccntroids of the d and / spin- 
orbit doublets is larger in the NL3 parametrization as 
compared with the NLl one by approximately 400 keV. If 
one corrects the NLl results by this energy gap, one gets 
the results indicated as NLl'^'"' in Fig. 7. The NLl™'' 
results are very close to the NL3 ones, which strongly 
suggests that the difference between the NLl and NL3 
results is predominantly due to different relative energies 
of the / = 3 and I = 2 states in these parametrizations of 
the RMF Lagrangian. 

Ref. [10] has attributed the fact that the NLl and NL3 
parametrizations differ in the description of the absolute 
AE value (Fig. 7) to the magnitude of the iso-scalar spin- 
orbit potential. The current investigation does not sup- 
port this interpretation. 

These results suggest that instead of readjusting the 
isoscalar strength of the spin-orbit interaction as it was 
done in Ref. [1], one can attempt to readjust the coupling 
constants of the DFT terms influencing the relative en- 
ergies of the I = 2 and I = 3 states with the same effect 
on AE. Indeed, 5% reduction of the isoscalar strength of 
the spin-orbit interaction introduced in Ref. [1] reduces 
AE by ^ 350 keV and this change in AE almost does not 
depend on nucleus (Fig. in Ref. [39]). On the other hand, 
the CRMF results suggests that the same eff'ect can be 
achieved if the relative distance of the d and / states is 
modified. Indeed, the AE values obtained in the CRMF 
calculations decrease by approximately the same amount 
on going from the NL3 to NLl parametrization of the 
RMF Lagrangian and this decrease only weakly depends 
on the nucleus (Fig. 7). 



3. The n-dependence of {AE*^ - AE^'"')norm 

The {AE*^ — AE'^^P)norm quantity shows pronounced 
dependence on n (Fig. 8a) and its trend (if normalized 
to a single nucleus) almost does not depend on the RMF 
parametrization. In order to understand which ingre- 
dients of AE^*-^ contribute into this n-dependence, the 
variations 5Ei = AEi{nucleus) — AEi{^'^V) of different 
(i-th) terms contributing to AE^^ are studied below. 
Contrary to Fig. 8a, ^''V is selected as a reference in or- 
der to get a picture less disturbed by large fluctuations 
of some variations in the vicinity of ^^Ca (Fig. 11). The 
6Eto = S{ETo{d^l2 f-ii-i ) - ETo{fy2}) variation is ob- 
tained in the deformed CRMF calculations (from Fig. 
8b), while other variations shown in Fig. 11 are calcu- 
lated in spherical RMF calculations. Thus, I effectively 
employ the approximation given in Eq. (2) assuming that 
the deformation polarization effects are the same both in 
theory and experiment. All the results presented here are 
based on the calculations with the NLl parametrization, 
but it was checked that the NL3 results are similar. 

The S{{AE*'' - AE''''P)norm) variation indicates that 
the difference between the calculated and experimental 
AE values decreases with decreasing n. Note that for 
a given n it is almost constant indicating only weak 



isospin dependence of this variation. The largest changes 
as a function of nucleus amongst the calculated varia- 
tions are seen in the energy gap between the centroids 
of the d and / spin-orbit doublets (the curve denoted as 
"SE{1 — centroids)" in Fig. 11). It has the same trend 
as S{{AE*'^ - AE'^^P)norm) as a function of n. For a 
given n, it shows very large dependence on isospin. The 
second largest variation is seen in the spin-orbit split- 
ting of the /7/2 — /5/2 spin-orbit doublet (the curve de- 
noted as ''5Ei,{fT/2 - /5/2)/2" in Fig. 11). The factor 
1/2 is used in this variation since only one half of the 
total variation of spin-orbit splitting contributes into the 
h/2 - splitting (see Sect. IVC 1). The 5Eis{fr/2 - 
/5/2)/2 variation has the wrong trend as compared with 
S{{AE*'^ - AE^^P)„)- The SE^d^/^ - 4/2)72 varia- 
tion of the spin-orbit splitting in the d^/2 ~ doublet 
is quite small. It has the correct trend as compared with 
5{{AE'h _ AE^='P)„). 

The SE{fY/2—d3/2) variation of the /7/2— ^3/2 splitting 
approximately satisfies the relation 

5E{f'ji2 — £^3/2) = 5E{1 — centroids) -|- 

5Ei,{fy2 - /5/2)/2 + 5Eu{d^,2 - d3/2)/2 (7) 

The isospin dependencies seen in 5E{1 — centroids) and 
5Eis{fT/2 — /5/2)/2 act is opposite directions, thus, re- 
ducing the isospin dependence of 5E{f-j/2 — ^3/2) as com- 
pared with the one of 5E{1 — centroids). However, the 
SE{f-j/2—d3/2) variation (Sect. IVC 1) cannot completely 
account neither for absolute value nor for isospin depen- 
dence (for a given n) of the S{{AE*'^ — AE'^^P)norm) vari- 
ation. 

Only when the 5E{f'ji2 — ^3/2) variation is combined 
with the 5Eto variation due to NM by 

SE^um ^ 5E{fy2 - ds/2) + 6EtO (8) 

a better description of the 6{{AE*'^ — AE'^^P)norm) vari- 
ation emerges. For a given n, the isospin dependence 
of the SiiAE^f" - AE^^P)„o™) is well described by 
(5^;""™. The absolute value of d{{AE^'' - A£:'=^P)„orm) 
for the Ti nuclei is well described by dE^""^. How- 
ever, for the Ca and Sc nuclei, the difference between 
these two quantities reaches 30% of the absolute value of 
5{{AE*'' - AE'''=P)norm)- The part of this discrepancy is 
definitely related to the limitations of the approximation 
given by Eq. (2). 

Thus, the current study clearly shows that the modi- 
fications of the relative placement of the states with dif- 
ferent angular momentum I, the spin-orbit splittings and 
time-odd mean fields on going from '^'^Y to ""^^Ca con- 
tribute into the n-dependence of the difference between 
the calculated and experimental AE values (the AE*^ — 
AE^'^P)norm quantity). Previously, this n-dependence of 
{AE*'^ — AE'^^Pjnorrm cxprcsscd in a different form (Fig. 
1 in Ref. [10], has been solely attributed to the deficiency 
of the iso- vector term of the spin-orbit interaction [10], 
but the current investigation does not support such an 
interpretation. 
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FIG. 11: (Color online) The variations SEi of different terms 
contributing to the S{AE*''' — AE'^'''^)„orm.) variation, see text 
for detail. 



V. THE ENERGY SCALE AND THE 
EFFECTIVE MASS OF THE NUCLEON 

The terminating states are expected to be of predom- 
inantly single-particle nature [14, 23]: the (^3/2 /"/2^ ter- 
minating states are obtained from the terminating 
state by particle-hole (p-h) excitation from the 1^3/2 state 
into the /y /2 state. The energy of this p-h excitation de- 
pends on the energies of above mentioned states, and, 
thus, it is affected by the energy scale of the single- 
particle spectra which is related to the effective mass 
m*{kF)/in of the nuclcon at the Fermi surface. 

In the RMF theory, the spin-orbit interaction is effec- 
tively scaled by the effective mass of the nucleon (Ref. 
[40]), and that is a reason why experimental data on 
spin-orbit splittings are well described in the calcula- 
tions [5, 36]. This scaling also explains why the spin- 
orbit splittings of the IP3/2 — lpi/21 l'^5/2 ~ l'^3/2 £ind 
1/7/2 ~ 1/5/2 spin-orbit partner orbitals are almost the 
same in the RMF and the Nilsson potential calculations 
(Fig. 5). Note, that the Nilsson potential is characterized 
by the effective mass m*(fci?)/m ~ 1 which is typical for 
experimental density of the quasiparticle states. Only 
in the case of the 2^3/2 — '2pi/2 doublet, the spin-orbit 
splitting is smaller in the RMF calculations. 

On the other hand, the energies of the centroids of the 
spin-orbit doublets are stretched out in the RMF cal- 
culations as compared with the Nilsson potential: the 
difference between the RMF and Nilsson centroid ener- 
gies increases on going away from the Fermi level (Fig. 5). 
Thus, the stretching out of the single-particle spectra due 
to low effective mass of the nucleon shows up mostly for 
orbital motion of particles and affects the relative place- 
ment of the levels with different angular momentum I. 
The origin of this problem has been discussed in Sect. 
IV C 2. 



The effective mass of nucleon at the Fermi surface 
(Lorentz mass in the notation of Ref. [41] for the case 
of the RMF theory) is m*{kF)/m - 0.65 in the RMF 
theory [5], ~ 0.7 in the case of the Hartree-Fock (HF) 
approach based on the Gogny forces [2], and varies in 
the range 0.6 — 1.0 in the HF approach based on the 
Skyrme forces [2] revealing much larger flexibility of this 
type of the DFT with respect of effective mass. As a con- 
sequence of low effective mass, the calculated spectra are 
less dense than the experimental ones: the well known 
fact in non-relativistic and relativistic models both for 
spherical [20, 36, 42] and deformed systems [43, 44]. This 
study shows that the AE^"-^ quantity differs from the 
fj/2 — energy gap in the spherical single-particle 
spectra only by the effects of time-odd mean fields and 
deformation polarization effects (Sect. IIIB2). These 
facts pose an open problem on how to compare the ex- 
perimental data on AE with the results of the DFT cal- 
culations (especially, those with low effective mass) since 
the experimental data on AE"^^^ is expected to be char- 
acterized by m*{kp)/m ~ 1. The implicit assumption 
used in Refs. [1, 10] that the DFT reproduces the em- 
pirical AE values relatively well, say within ^ 10% [39], 
may be too optimistic especially for the DFT with low 
effective mass. 



VI. CONCLUSIONS 

In conclusion, the following results were obtained in 
the study of band termination within the DFT frame- 
work: 

• At band termination, the NM does not modify 
neither total angular momentum nor the expecta- 
tion values of the single-particle angular momenta 
< jx >i of the single-particle orbitals. NM pro- 
vides an additional binding to the energies of the 
specific configuration and this additional binding 
increases with spin and has its maximum exactly 
at the terminating state. This suggests that the 
terminating states can be an interesting probe of 
the time-odd mean fields related to NM provided 
that other effects can he reliably isolated. 

• The realization of the TS-method in Refs. [1, 10] is 
based on the analogy with simple form of the Nils- 
son potential which allows to neglect the deficicnccs 
in the relative placement of the states with differ- 
ent angular momentum I. This approximation is 
not valid for terminating states in the A ^ 44 mass 
region in modern and most frequently used versions 
of the Nilsson potential 

• The impact of the relative placement of the states 
with different angular momentum I on AE^'^ is 
also clearly visible in the RMF calculations. The 
difference in absolute AE^'-' values obtained in 
the CRMF calculations with the NLl and NL3 
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paramctrizations of the Lagrangian is defined by 
the different relative energies of the I = 3 and I — 2 
states in these parametrizations. The modifications 
of the relative distance of the states with different 
angular momentum / on going from "^^V to ^^Ca 
contribute into the n-dependence of the difference 
between the calculated and experimental AE val- 
ues (the AE*'^ — AE'^^P)norm quantity) in addition 
to the ones due to the spin-orbit interaction and 
time-odd mean fields. 

The detailed analysis of the TS-method in the RMF 
framework reveals the picture which is more complicated 
than the one suggested in Refs. [1, 10]. The relative 
placement of the states with different angular momen- 
tum /, defined by the properties of central potential, has 
to be taken into account in addition to the DFT terms 
discussed in these references when the AE quantity is 
analyzed. Considering the similarities of the RMF the- 
ory and SDFT, it is very likely that these conclusions are 
also valid in the SDFT framework. The current inves- 
tigation calls for a detailed study of the impact of the 
relative placement of the states with different orbital an- 
gular momentum I on the AE^^ quantity in the SDFT 
framework. 

Existing results for superdeformed bands in '^^S [45, 46] 
and low-spin states in odd mass nuclei [32] point to the 



time-odd mean fields as a major point of the difference 
between SDFT and RMF. For example, the additional 
binding due to time-odd mean fields and the energy sep- 
aration between different signatures of the SD bands are 
considerable stronger in SDFT as compared with RMF 
[45]. The current study clearly shows that the corre- 
lations induced by time-odd mean fields are large: ad- 
ditional binding due to NM reaches 4 MeV for n = 3 
(Fig. 9), which is by order of magnitude larger than those 
seen before in the RMF calculations at low spin. It also 
has a considerable impact on AE^'^: AEto ~ 1-2 MeV 
(Fig. 8b). These results call for a comparative study 
of time-odd mean fields in the Skyrme DFT and RMF 
frameworks. Such study is necessary in order to make 
a significant progress towards a better understanding of 
the role of time-odd mean fields. The work in this direc- 
tion is in progress and the results will be presented in a 
forthcoming manuscript [32]. 
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